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!• GENERAL EQUATION OF THE VIBRATING 

MEMBRANE. 

The problem of the vibrating membrane has been attacked 
by previous writers from the standpoint of the shape of the 
bomidaries of the membrane, rectangular, circular, parabolic, 
eUiptic. No attempt has been made to derive a single equation 
which contains all of these. Here we shall derive such an equa- 
tion by the use of a iteneral system of curvilinear codrdinates. 
These codrdinates shall be restricted only by the condition that 
the resulting equation shall have a solution in the form of a product 
of functions of the independent variables involved, each function 
depending only upon a single variable. 

If a stretched elastic membrane is fastened at the edges and is 
made to vibrate at right angles to the plane of the membrane, 
its equation of motion in rectangular codrdinates is, 

where a is a physical constant depending upon the mass per 
imit area and the tension along any line in the surface of the 
membrane. In (1) put V ^ T-u, where 7 is a function of i only 
and u a function of x and y only. The equation then separates 
into the two parts, 

(2) ^ + Ip+^ = 0, and ^+a'K*T^O, 

where K is any constant. Solutions of the latter are 7 = sin aKi 
and r = cos aKt. 

Now suppose X and y are determined as functions of two new 
unrestricted variables, x = /i(a, p), y = f%{a, p) and express 
dht/da? + dhi/dy^ in terms of these new variables, 

d3^'^dy^^d^l\di) '^\dy) y dP^[\dx ) 
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(3) +m]+^\p.f+p.f'\ 

^ . \oy / J dadp\_dx dx dy ay J 



Assume now that 
(4) 



x+iy^f(a + ifi), 
X — iy = /(a — iff). 

This assumption will greatly simplify (3). Putting this simpli- 
fied result in (2) we obtain, 

By means of the orthogonal relation 

(g)*+(^)'=-[(ii)*^(rj] 

this equation becomes 

c) ^+P=-«-[(l)'+(E)"]- 

The right member of this may be written in a more convenient 
form. From (4) we may derive the relation, 

which when put in (6) gives, 

(8) ^ + Qp2^-IPu'f\a + ip)'f'{a-'ifi). 

On the assumption that this equation has a solution of the form 
u s= 6{a) 'H(p), it is necessary that the relation (9) should hold. 

(9) / '{a + ip) •/ \a -iP)^ E{a) + F(fi) 

where E and F are functions which we shall determine now. 
By differentiating first with respect to a and then with respect 
to p we obtain the equation 

(10) /'"(a + ip) 'f\a -ip)- S'"(a - ip) -/'(a + tf ) « 0. 
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This equation has the solutions 

(11) f\a + iP) = BsinAia + iP) and C cos A{a + ip), 

A, B, C are constants. 

These are the functions then which satisfy the assumption in (9). 
Substituting in (9) we obtain^ 

E{a)+ F(fi) = B^ sin A{a + tjS) -sin A(a - i/3) 

= B\siD? Aa + sinh« Aff) 
and (8) becomes, 

g + P + IPBMw^^ AP + sin* Aa) = 0. 
Hence (1) may now be written, 

(12) [B«(3inh« A^ + sin* Aa)V^ "^ "Aa^ + ajS^ j ' 

That the assumption in (9) really has resulted in changing our 
curvilinear coordinate system to families of orthogonal confocal 
ellipses and hyperbolas may be seen from (11) and (4). And the 
method here employed shows that (12) is the most general 
equation obtainable for the vibrating membrane which has for 
its solution the function V = T{()*G{a)'H(fi). 

2. DERIVATION AND DISCUSSION OF THE GENERAL 

EQUATIONS OF THE CONFOCAL CONICS 

USED AS COORDINATE SYSTEMS. 

From the equations of the families of the ciu'ves of the orthog- 
onal ciu'vilinear coordinate system in which (12) is expressed, 
we purpose now to derive the following systems as special cases^ 

(a) System of rectilinear lines. 

(6) System of concentric circles and radial lines. 

(c) System of confocal parabolas. 
In order to do this, however, care must be exercised in selecting a 
general solution for (10). Take this solution in the form 

f{a+ip) ^ BsinA{a + iP) - C cos A{a + ifi). 
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The system of rectilinear lines obtained in (a) may jtlso be 
derived from^^the system of confocal parabolas. Put A^ = p 
in (22) and Ao? = 9 in (23) and change the origin to the point 
(o, 0). 
(24) f^2p{x+a) + p', 

(26) j» = - 2q{^ + a) + if. 

Now if we let 2pa = m* in (24) and then let a = 00 in the 
negative direction while m remains finite, we obtain 

(26) y = ± m, 

and this condition can be brought about by having /3 ^ as 
a =s 00 in such a way that 2A^'a = m^ for A is finite. 

If we let a — 9/2 = — n in (26) and then let a = 00 while n 
remains finite, 

(27) x^n. 

And this condition can be brought about by having a « 00 as 
a = 00 in such a way that {Ao?l2) — a always reniains finite. 
For those parabolas of the family in (23) in which {Ac^/2) > a 
the parameter n is positive and for those in which {Ac^/2) < a 
the parameter n will be negative. 

If we had changed the origin of the curves in (14) and (16) to 
the right focus and then sent the left to infinity in the negative 
direction, we would have obtained 

f = Acf{Ao^ + 2x) and f = A0^(A0^ - 2x) 

as the equations of the families of parabolas. From these the 
families of rectilinear lines may be obtained by changing the 
origin to the point (— a,0) and then sending the focus to infinity 
in the positive direction. 

3. DISCUSSION OF THE GENERAL DIFFERENTIAL 
EQUATION OF THE VIBRATING MEMBRANE. 

The same values of the constants used for the derivation of 
the three special systems of curvilinear coordinates from the 
general system should, when put in the general equation of the 
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vibrating membrane, with the time element removed, give three 
special equations expressed in the above mentioned coordinates. 
With the aid of (13) equation (6) may be written 

(28) ^ + ^ + li?C*i4^[cosh« (^l8 + D) 

- cos* {Aa + J?)] = 0. 
This separates into, 

(29) ^H- [K^CA^ cosh« {Afi + D) - iP\Y = 0, 

(30) ^ - [K^C'A^ cos* {Aa + B) - M']X = 0. 

The family of horizontal lines, y =» =b Fj9, were obtained from 
the family of confocal ellipses by having C = oo and ^4 » in 
such a manner that C-A — F and Z> = 0. If these values are 
put in (29), it becomes 

(31) ^+(ii?-3f«)y=0 where M^j. 

The family of vertical lines « = =F Fa were obtained from the 
family of confocal hyperbolas by having C = oo and ^ » in 
such a way that C'A ^ F and B = ir/2. If these values are 
put in (30) it becomes 

d^X — — M 

(32) ^+3f«Z=0 where M^j. 

And these are the equations into which 

separates. 

The family of concentric circles, a:* + j* = F^^^, were ob- 
tained from the family of confocal ellipses by having Z) = oo 
and C ^ in such a way that C cosh D = F and C sinh D ^ F. 
Moreover 

d^Y <py dY 
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in this instance. If these values are put in (29) it becomes 

K^KF and M = ^- 

The family of radial lines, y ^ ±x tan Aa^ were obtained 
from the family of confocal hyperbolas by having C » and 
B = 0. If these together with 

ip = tan-^-= Aa 

X 

are put in (29) it becomes 

(34) ^ + Jf2Z = where J^ = J- 

And these are the equations into which the equation which de- 
fines the function u in circular coordinates^ 

Shi , 1 dti , 1 3*ti . -^ 

separates. 

It will simplify matters somewhat, in obtaining the equations 
in parabolic coordinates, if we adopt the equation 

(35) ^+J^+ IPCAM^oh* (Afi + D) + sin* (Aa + B)] = 
in place of (28). This separates into 

(36) ^+ [E?C^A^ sinh« (Ap + D) - AP\Y = 0, 

(37) ^ + [IPC'A^ sin* {Aa + B) + AP\X = 0. 

The family of confocal parabolas, y^ = A^{A^ + 2a;), if ere 
obtained from the family of confocal ellipses by having A ^ 
and C = 00 in such a way that C'A^ ^ A and Z> = 0. If these 
values are put in (36) it becomes 

(38) ^+{l^e^-iP)Y^Q, where 1f-^K}^AK 



THE VIBRATING MEMBRANE. 9 

The family of confocal parabolas, y^ = Ao?{Ao? — 2x) were 
obtained from the family of confocal hyperbolas by using the 
same values of the constants as above and £ =» 0. If these are 
put in (37) it becomes 

(39) ^+(ifc»a*+iP)X = where *« = ii?.I«. 



4. SOLUTION OF THE EQUATION IN RECTANGULAR 

AND CIRCULAR COORDINATES. 

The complete solution of the equation of a stretched elastic 
rectangular vibrating membrane is obtained by an extension of 
Fourier's Theorem and may be found in the works of any of 
the writers on Partial Differential Equations. The complete 
solution representing the vibrations of a circular membrane is 
obtained by the use of Bessel's Functions. If the membrane 
vibrates so that it has circular symmetry about the axis per- 
pendicular to the plane of the boundary at the origin^ the Jo 
functions furnish the solution. If the mode of vibration is a 
function of both the coordinates r and ip then the Jn functions 
furnish the desired solution. And a discussion of the problem 
may be found in any of the works on Bessel's Functions. What 
we believe, however, to be a contribution to the literature of 
Bessel's Functions, is a method of obtaining the following funda- 
mental relations between these functions: 

dJn{x) n 

(a) —^ = - Jnix) - J^lix), 

(b) 2 ^^ - J^iix) - J^^ix), 

(C) ^^^J^l(^)-ljn{xh 

(d) ^Jn{x) = J^iix) + J^iix), 
X 



ie) 



da? 



= (^^-l)va!/.(«). 
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The usual method of obtaining these relations is to consider the 
functions entirely independent of the equation of which they 
are solutions. Since the functions are defined by a differential 
equation it was thought that the equation itself ought to yield 
these fundamental relations. 

Starting with Jn{x) as a solution of Bessel's Equation, 

and substituting z = ^v^ we obtain v = x'^Jnix) as a solution of 

d^v 2n+ 1 dv 

Substituting this value or « in (41) 

cP 2n + 1 d , 

(42) ^ (pr*J„x) -\ — ^ {?r*Jnx) + ar» J,ie = 0. 

The corresponding equation for J»\.iX is 

(43) d;? (* -^-H-i*) + -T- ^ (^ '•^■H-i*) 

Differentiate (42) with respect to x, add to (43) multiplied by x 
and arrange terms, 

A solution for this is 

d 
(45) ^ (ar*/nir) + x-^Jn^ix - 0. 

A corresponding relation between JnX and J»-ia; may be 
obtained. Put z = x'^v in (40) 
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d?t 2n — 1 dt^ 

A solution for this is v — 3^JnX. Put this in (46), 

(47) ^^_?^'^^+^/^.0. 

The corresponding equation for Jn^ix is 

(48) ^(»^»-i^) _ 2!^ j(£r^ ^ ^_,^^^^ ^ Q 

Differentiate (47) with respect to x, subtract (48) multiplied by 
X from it and arrange terms 



^ [^ («-/,») -a^J^ia; J 



(49) 



2n— Idfi,,, , 1 



A solution for thb is 

d 
(60) ^ (afJnx) - a*J,^ia; = 0. 

Relations (a), (6), (c), ((f) may be obtained from (45) and (50). 
Relation (e) may be obtained from (40) by the substitution, 
s = t/-^x, where » = -^IxJ^x. 

5. SOLUTION OF THE EQUATION IN PARABOLIC 

COORDINATES. 

The equation of the vibrating membrane written in parabolic 
codrdinates is 

(51) («* + ^a^=«'id^ + d^j- 

This separates into 
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and 

(52) ±-^+(j^c? + \)X = 0, 

(53) ^+(*^i8»-X)y = 0. 

Dr. H. Weber in his Partiellen Differentialgleichungen^ Vol. 2, 
page 256, finds two imaginary solutions for (52) as limiting forms 
of the hypergeometric series. He also finds two real solutions 
in definite integral form. Several new solutions for special 
values of X will now be obtained. 

Repeated differentiation of (52) will point the way to the 
construction of equations which when differentiated 1, 2, 3 • • - n 
times will give (52). The nth equation of this series can be 
solved if X has the special value X = iki^n + 1), where n is a 
positive integer. And from this solution a solution of (52) may 
be obtained. It is 

(54) -X"" C(V «)«•'. ^iC*-^*). 

From this a series of functions depending upon n may be ob- 
tained. For n a we have the very simple one, 

(55) Z= C(V «)-**•'. 
Weber's two solutions f or X = tJk become 

Zi = C,(Vc)*^' • r 1 + (- »*a*) + ^ (- iko?)* 

(56) •- /'^ 

+ r2:3^- *«*>•+ -J' 

X, = C,( V e)'*^V (- i*a») [ 1 + j^ (- tfco*) 
(57) 

Equation (56) reduces to Zi = Ci( V 6)^**** and this agrees with 
(55). The series in (57) may be summed by means of the 
integral 
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j^ g-'i^l xr^dx = 2e-*Vx[l + j^« + 1^:5 «* +•••], 

which is obtained by repeated integration by parts of the left 
member. Applying this to (57) we obtain, 



None of the above results lend themselves to the solution of the 
problem of the vibrating membrane because they are imaginary. 
For X = 0, however, two real solutions of (52) are 



(59) 



Xi s an I 1 — ^-^ A -— ^^ 1- • • • I 

Ai oo|i 3-4^3-7-4-8 3-711-4-812 ^ J' 



Z, = o,(i/^)ri-^V 



iV 



(60) •- jfc»a» 

4-812.5-913"^ 



] 



The corresponding F-functions for (5S) may be found from these 
by replacing a with /3. These functions then may be used in 
solving the problem of the vibrating membrane with parabolic 
boundaries when the mode of vibration is such that X »=: 0. 

The first two roots of Zi = and Zj = in (59) and (60) 
have been computed by the writer. They are 

#^ r nl^i=*-013--., . ^ .fxi= 5.563..., 
forZx=0J^^^j0.246..., '^^^« = M ^2 = 11.818.... 

The roots for Fi = and 72 = are identical with the corre- 
sponding roots of Zi = and X% = 0. These roots are the 
values of the parameter which give the nodal lines for this par- 
ticular mode of vibration. 

Weber^s two imaginary solutions for X = are, 

Xi = (v«)'*-ri + i(- 1*«^ + 4^(- *«•)• 

(61) 
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Xt = ( Ve)**-' V (- ifto*) I 1 + |(- iJco?) 
(62) 7 11 1 

+ 5T2-I <- *«•)' + 5:2^1 (- ^'- • \ ' 

Equation (61) reduces to (59), a real solution, if oq = 1. Equa- 
tion (62) reduces to (60), a real solution, if ai = V— t. These 
results may be obtained by performing the indicated multipli- 
cations in each case. 



6. REVIEW OF THE WORK PREVIOUSLY DONE ON 

THE SOLUTION OF THE EQUATION IN 

PARABOLIC AND ELLIPTIC 

COORDINATES. 

Parabolic coordinates: The real solutions of the equation 

obtained by Dr. H. Weber* and referred to previously are 

(63) -STi = jT *-*(! - s)-^ cos [ft^a*(* - i) + ^log j^-^]** 

(64) Z, - jr^-*(I ~*)-*cos[A2a*(*-i) +^logj-£:^]A- 

The corresponding functions of j3 may be obtained from these by 
changing a into j3 and X into ^ X. Using these solutions Weberf 
discusses the vibration of regions bounded by parabolic nodal 
lines. He proves the theorem, " If a function of u different from 
zero, defined by the equation 

d^u Shi 

is continuous and its first differential quotient also, inside a 
bounded region but vanishes on the boundary, then k^ must be 
real and positive." Such regions are of three kinds, 

* ''Partielle Differentialgleichungen/' Vol. 2, page 256. 
t MathemaiUche Anruden, Vol. 1, 1868. 
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(a) Bounded by two parabolas one from each family. 

(6) Bounded by two parabolas from one family and one from 
the other. 

(c) Bounded by two parabolas from one family and two from 
the other. 
For (a) the boundary conditions are Z = for a = =t ai and 
y = for j8 = =t /3i. Under these conditions four particular 
solutions may be found, u^ Xi'Yi, w = Zj • 72, u^ Xi- Yt, 
K = Z2 • Yi. The last two however will not satisfy the condition 
that u shall be single valued within the region. Two separate 
systems of solutions then may be built from li = XfYi and 
u ^ Xf Y%. In each solution the constants X and k appear and 
may be determined as the roots of the two transcendental equa- 
tions involved. 

For (6) the boundary conditions are X = for a = ai, a = aj 
and 7=0 for jS = ± j8i. Two types of solutions will fit 
these conditions, u = Yi{AiXi + AiX%), u = Yi(AiXi + AtXi). 
These give rise to two systems of solutions. The constants in 
each solution are A1/A2, k, X and they may be determined from 
the three transcendental equations involved. 

For (c) the boundary conditions Are Z = for a = ai, a = a% 
and y = for jS = j8i, jS == ft. One system of solutions only 
can be used here, u = {AiXi + A2X2)(BiYi + BtYt) and there 
are four transcendental equations to determine the constants 
A1IA2, B1/B2, X, *. 

M. Hartenstein* obtains solutions for the equation 

dhi , Shi , -, 

in parabolic coordinates for negative values of k^ by trans- 
forming the solutions of the equation in circular coordinates 
(Bessel's Functions) into their corresponding functions expressed 
in parabolic coordinates. This is a part of a general discussion 
of the equation expressed in the four systems of coordinates used 
in the third part of this paper. The results obtained contribute 
nothing to the discussion of the problem of the vibrating mem- 
brane. 

* Archw der MathemaUk und Phynk (2), Vol. 14. 
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Elliptic coordinates: The equation of motion of the vibrating 
membrane expressed in eUiptic coordinates is^ 

(65) (co3h»^-cos^a)^=a«(^+a^j 

and separateis into 

df'T 

and 

(66) dc? "^ dip "^ *^(cosh* /3 - cos* ct)u = 0. 

This last equation separates into^ 

(67) ^-(i*cos2a-X)Z=0, 

^^^^ ^^ ^^ ^^^ jS - X)y = 0. 

Dr. Heine* shows that (66) is a limitins: form of the differential 
equation of the Lam6 Functions in the same way as the differ- 
ential equation of the Bessel's Function is a limiting form of 
the differential equation of the Spherical Functions. He solves 
(67) by assuming that the constant X depends in a definite way 
(as a root of a known transcendental equation) on the constant k 
and the solutions periodic functions of a. Four classes of func- 
tions are found to be solutions of (67) under these assumptions. 
They are, 

(a) Series of cosines in ascending even multiples of a. 

(b) Series of cosines in ascending odd multiples of a. 

(c) Series of sines in ascending odd multiples of a. 
(cO Series of sines in ascending even multiples of a. 

E. Mathieuf in discussing the vibratory movement of an 
elliptic membrane, uses periodic functions as solutions also of 
the above equations. He shows that the general solution of 
such linear differential equations of the second order as (67) 
and (68) is composed of two parts, one zero for the zero value of 
the argument and the other a maximum. If the constant X has 

* "Handbuch der Kuglefunctionen/' page 401. 
t Jwmak de LiaumUe (2), Vol. 13, 1886. 
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the form X » y* + 6i* + ci* + • • • where 9 is an integer and 
b, e, etc. constants determined by the periodic property, the 
solutions of (67) and (68) may be written as periodic functions 
of a and /3 in ascending powers of 1^, and they become zero g 
times from to t. For thespecial case when ik » these solu- 
tions become Xi » sin (ga), JTs » cos (ga), Yi » sin (jgifi), 
Yt » cos, igiP)* Solutions for these equations are also given in 
series according to the ascending powers of sin and cos of a 
and /3. In this case the general solution of (67) is the sum of 
solutions one even and the other odd in sin a and one even and 
the other odd in cos a. These results agree with those given 
by Heine. From these solutions two types of vibration are 
possible. And some simple cases of hyperbolic nodal lines are, 
the major axis alone, both together, two asymptotes of the same 
hyperbola, two asymptotes of the same hyperbola with the major 
and minor axes. The elliptic nodal lines are obtained from 
7 1 = and Fs » as foQows: If fi ^ B on ihe boundary of 
the memb]:ane when 7 « 0, the constant X may be found from 
Y{B, X) ^ 0. If Xi, Xs, Xg, • • • are the roots of this equation in 
order of increasing magnitude, Y(fi, Xn) « will give, through 
its roots in /3, the parameters of the elliptic nodal lines. This 
equation is shown to have n — 1 roots less than B, and there- 
fore the number of elliptic nodal lines is n ^ 1. 

F. Pockels* has investigated the function defined by (66) for 
the region bounded by any two ellipses, P ^ Pi and P = fit and 
any two hyperbolas, a « ax and a » as under the supposition 
that u = on all four sides. He proves the foQowing theorem: 
''There is for such a region a doubly infinite series of normal 
functions which satisfy the boundary condition u «= on all 
four sides, and of which a definite number (m ~ 1) of elliptic 
and a definite number (n — 1) of hyperbolic nodal lines are 
characteristic, (m = 1, 2, • • •, oo, n » 1, 2, • • •, oo)." By cut- 
ting the bounded region along the line connecting the foci, and 
assuming that u and its first derivative are continuous along this 
new boundary, the new region thus formed may be regarded as a 
simply connected Riemann surface. This makes possible the 
consideration of regions bounded by less than four curves, as 

* "IJeber die Partielle Differential^eichttng An + )b«tt « 0." 



